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Recently, several white dwarfs have been proposed with masses significantly above the Chan¬ 
drasekhar limit, known as Super-Chandrasekhar White Dwarfs, to account for the overluminous 
Type la supernovae. In the present work, Equation of State of a completely degenerate relativis¬ 
tic electron gas in magnetic field based on Landau quantization of charged particles in a magnetic 
field is developed. The mass-radius relations for magnetized White Dwarfs are obtained by solving 
the Tolman-Oppenheimer-Volkoff equations. The effects of the magnetic energy density and pres¬ 
sure contributed by a density-dependent magnetic field are treated properly to find the stability 
configurations of realistic magnetic White Dwarf stars. 
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I. INTRODUCTION 

Ultrahigh magnetic fields in nature are known to be 
associated with compact astrophysical objects namely 
white dwarfs, neutron stars and black holes. Of these, 
the largest magnetic fields are found on the surfaces of 
magnetars, Anomalous X-ray Pulsars (AXPs) and Soft 
Gamma Repeaters (SGRs), certain classes of neutron 
stars, with an order of magnitude of 10 15 gauss. Recently, 
a strong magnetic held of the same order of magnitude 
as that of a magnetar has been found at the jet base of a 
supermassive black hole PKS 1830-211 pj. These strong 
magnetic fields drastically modify the Equation of State 
(EoS) of a compact star and its stability. Hence, study¬ 
ing the EoS and equilibria of compact stars in presence of 
high magnetic fields is an important and rapidly growing 
field of research in theoretical astrophysics. Magnitudes 
of magnetic fields of white dwarfs are constrained by the 
virial theorem: 


The recently observed peculiar Type la supernovae, 
e.g. SN2006gz, SN2007if, SN2009dc, SN2003fg, [2^3 
with exceptionally high luminosities do not fit with the 
explosion of a Chandrasekhar mass white dwarf. More¬ 
over, it has been seen that there is a correlation between 
the surface magnetic field and the mass of white dwarfs. 
The magnetic white dwarfs seem to be more massive than 
their nonmagnetic counterparts []}. Lastly, predictions 
from the luminosities reveal that the progenitor white 
dwarfs had masses significantly higher than the Chan¬ 
drasekhar limit. It seems that the Chandrasekhar limit 
may be violated by highly magnetized white dwarfs. To 
account for these facts, we have calculated theoretically 
the masses of white dwarfs in presence of such high mag¬ 
netic fields in the general relativistic formalism. 

II. COMPLETELY DEGENERATE IDEAL 
FERMI GAS AND EOS FOR NON-MAGNETIC 
WHITE DWARFS 
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Here, B , Bq, M, M 0 , R , Rq are the magnetic field, 
mass and radius of the white dwarf and sun respectively. 
Using Bq = 2 x 10 s gauss, M = 1.4 Mq and R = 0.0086 
Rq, we get the order of magnitude as B max ~ 10 12 gauss. 
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We consider a relativistic, completely degenerate Fermi 
gas at zero temperature and neglect any form of interac¬ 
tions between the fermions. By the Pauli exclusion prin¬ 
ciple, no quantum state can be occupied by more than 
one fermion with an identical set of quantum numbers. 
Thus a noninteracting Fermi gas, unlike a Bose gas, is 
prohibited from condensing into a Bose-Einstein conden¬ 
sate. The total energy of the Fermi gas at absolute zero is 
larger than the sum of the single-particle ground states 
because the Pauli principle implies a degeneracy pres¬ 
sure that keeps fermions separated and moving. For this 
reason, the pressure of a Fermi gas is non-zero even at 
zero temperature, in contrast to that of a classical ideal 
gas. This so-called degeneracy pressure stabilizes a white 
dwarf (a Fermi gas of electrons) against the inward pull 
of gravity, which would ostensibly collapse the star into 
a Black Hole. However if a star is sufficiently massive 
to overcome the degeneracy pressure, it collapse into a 
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singularity due to gravity. While the pressure inside a 
white dwarf is entirely due to electrons, its mass comes 
mostly from the atomic nuclei. 


A. Completely degenerate free Fermi gas 

The non-interacting assembly of fermions at zero tem¬ 
perature exerts pressure because of kinetic energy from 
different states filled up to Fermi level. Since pressure 
is force per unit area which means rate of momentum 
transfer per unit area, it is given by 
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where m e is the rest mass, v is the velocity of the particles 
with momentum p and n p d 3 p is the number of particles 
per unit volume having momenta between p and p + dp. 
The factor 4 accounts for the fact that, on average, only 
ird of total particles n p d 3 p are moving in a particular 
direction. For fermions having spin A, degeneracy = 2, 
n p d 3 p = S7T ^ 3 dp and hence number density n e is given by 
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where pf is the Fermi momentum which is maximum 
momentum possible at zero temperature, Xf = is 
a dimensionless quantity and A e = is the Compton 
wavelength. The energy density e e is given by 


£ e = J En p d 3 p = J sj (p 2 c 2 + m 2 c 4 ) (5) 

which along with Eq.(3) turns out upon integration to be 
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B. EoS for non-magnetic White Dwarfs 

For the EoS for non-magnetic White Dwarfs, the pres¬ 
sure is provided by the relativistic electrons only and 
therefore, pressure P is given by 


P = Pe = ( 9 ) 

whereas for energy density £ both electrons (with its ki¬ 
netic energy) and atomic nuclei contribute, so that 

£ = £e+n e (m p +jm n )c = —ra-x{xF)+n e (m p +}m n )c 

''e 

( 10 ) 

where m n and m p are the masses of neutron and proton, 
respectively and / is the number of neutrons per electron. 
Commonly, electron-degenerate stars consist of helium, 
carbon, oxygen, etc., for which / = 1. To be precise, 
one should in fact also subtract n e (l + /) times binding 
energy per nucleon from the second term on the right 
hand side of the above equation. Obviously, this correc¬ 
tion is composition dependent and its contribution being 
quite small, e.g. in case of helium star it is about 0.7% 
to the second term, it is not considered in calculations. 
Since the kinetic energy of electrons in the above equa¬ 
tion contributes negligibly, the mass density for / = 1 
white dwarfs can be expressed in units of 2x10° gms/cc 
by multiplying number density of electrons n e expressed 
in units of fm -3 by the factor 1.6717305x 10 6 . 


III. LANDAU QUANTIZATION AND EOS FOR 
MAGNETIZED WHITE DWARFS 

Like the former case, here also we consider a com¬ 
pletely degenerate relativistic electron gas at zero tem¬ 
perature but embedded in a strong magnetic field. We do 
not consider any form of interactions with the electrons. 
Electrons, being charged particles, occupy Landau quan¬ 
tized states in a magnetic field. This changes the EoS, 
which, in turn, changes the pressure and energy density 
of the white dwarf. In addition to the matter energy 
density and pressure, the energy density and pressure 
due to magnetic field are also taken into account. It is 
the combined pressure and energy density of matter and 
magnetic field that determines the mass-radius relation 
of strongly magnetized white dwarfs. It should be em¬ 
phasized that protons also, being charged particles, are 
Landau quantized. But since the proton mass is ~ 2000 
times the electron mass their cyclotron energy is ~ 2000 
times smaller than that of the electron for the same mag¬ 
netic field, and hence we neglect it. 


A. Landau quantization and EoS for free electron 
gas in magnetic field 

In order to calculate the thermodynamic quantities like 
the energy density and pressure of an electron gas in a 
magnetic field, we need to know the density of states 
and the dispersion relation. The quantum mechanics 
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of a charged particle in a magnetic field is presented in 
many texts (e.g. Sokolov and Ternov (1968) Lan¬ 
dau and Lifshitz (1977) [§], Canuto and Ventura (1977) 
0 Meszaros (1992) [III). Here we summarize the basics 
needed for our later discussion. Let us first consider the 
motion of a charged particle (charge q and mass m e ) in a 
uniform magnetic field B assumed to be along the z-axis. 
In classical physics, the particle gyrates in a circular orbit 
with radius and angular frequency (cyclotron frequency) 
given by 


m e cv _l 
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where v± is the velocity perpendicular to the magnetic 
field. The hamiltonian of the system is given by 


H = 



( 12 ) 


where B = V x A with A being the electromagnetic vec¬ 
tor potential. To have magnetic field in z-direction with 
magnitude B one must have 


number k y . If periodic boundary condition is assumed 
k y can take values k y = where N is another integer 
and l x> l y , lz being the dimensions of the system. The al¬ 
lowed values of N are further restricted by the condition 
that the centre of the force of the oscillator Xq must phys¬ 
ically lie within the system, 0 < Xq < l x which implies 
0 < N < lxl y^£ ulc = qB f i x c l ' y . Hence for electrons with 
spin s and charge q = — |e|, the maximum number of 
particles per Landau level per unit area is . On 

solving Schrodinger’s equation for electrons with spin in 
an external magnetic field in z-direction which is uniform 
and static, Eq.(16) modifies to 

p 2 1 

E v,p z = vTiuic + ’ v = n -\- — + s z . ( 17 ) 

Clearly for the lowest Landau level (v = 0) the spin de¬ 
generacy g u = 1 (since only n = 0, s z = — ^ is allowed) 
and for all other higher Landau levels (u / 0), g u = 2 
(for s z = ±i). 

For extremely strong magnetic fields such that hui c > 
m e c 2 the motion perpendicular to the magnetic field still 
remains quantized but becomes relativistic. The solution 
of the Dirac equation in a constant magnetic field 0 is 
given by the energy eigenvalues 



and therefore 
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where the dimensionless magnetic field defined as Bn = 
B/B c is introduced with B c given by Tilu c = h = 

m e c 2 => B c = = 4.414 x 10 13 gauss. Obviously, 

the density of states in presence of magnetic field gets 

(14) modified to 


The operator p y commutes with this hamiltonian since 
the operator y is absent. Thus operator p y can be re¬ 
placed by its eigenvalue Tik y .Using cyclotron frequency 
uj c = -2^- one obtains 

° m e c 
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the first two terms of which is exactly the quantum 
harmonic oscillator with the minimum of the potential 
shifted in co-ordinate space by Xq = hk * . Noting that 
translating harmonic oscillator potential does not affect 
the energies, energy eigenvalues can be given by 
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where the sum is on all Landau levels v. At zero temper¬ 
ature the number density of electrons is given by 
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where PFiy) is the Fermi momentum in the zAh Landau 
level and v m is the upper limit of the Landau level sum¬ 
mation. The Fermi energy Ep of the uth Landau level is 
given by 


1 p 2 

E n,p z = (n+ -)hu>c + n = 0,l,2.... (16) 

The energy does not depend on the quantum number k y , 
so there will be degeneracies. Each set of wave functions 
with same value of n is called a Landau Level. Each 
Landau level is degenerate due to the second quantum 


E 2 f = p 2 F (v)c 2 +m 2 e c A (1 + 2vB d ) (21) 

and v m can be found from the condition [p F (i/)] 2 > 0 or 
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where e F = m F c z is the dimensionless Fermi energy and 
e Fmax = the dimensionless maximum Fermi en¬ 

ergy of a system for a given Bf> and v m . Obviously, very 
small Bo corresponds to large number of Landau levels 
leading to the familiar non-magnetic EoS. v m is taken to 
be the nearest lowest integer. Like the former case, if we 
define a dimensionless Fermi momentum xf(v) = PF — 
then Eqns.(20) and (21) may be written as 

Ue = (2nV\3 J2 ^XfM (23) 

' ’ e i /=0 

and 

£F=[4M + l+2 vB D y (24) 


or 


x F (y) = [e 2 p - (1 + 2 vB d )] ' 
The electron energy density is given by 
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where 


V’(^) = [ (1 +y 2 ) 1/2 dy = i[^ a/ 1 + z 2 + \n(z+y/l + z 2 )\ 
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The pressure of the electron gas is given by 


Pe 
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where 


rj(z) = z \/1 + z 1 —ij){z) = ^-[zy /1 + £ 2 —ln(z-f -\/1 + z 2 )]. 

(29) 


B. Magnetized White Dwarfs 

In the present case of magnetic White Dwarfs, the ex- 
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plicit contributions from the energy density eb = #- and 


pressure Pn = \eb arising due to magnetic field need to 
be added to the matter energy density and pressure as 
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IV. TOLMAN-OPPENHEIMER-VOLKOFF 
EQUATION AND THEORETICAL 
CALCULATIONS OF MASS-RADIUS RELATION 
FOR WHITE DWARFS 

If rapidly rotating compact stars were non- 
axisymmetric, they would emit gravitational waves 
in a very short time scale and settle down to axisym- 
metric configurations. Therefore, we need to solve 
for rotating and axisymmetric configurations in the 
framework of general relativity. For the matter and the 
spacetime the following assumptions are made. The 
matter distribution and the spacetime are axisymmetric, 
the matter and the spacetime are in a stationary state, 
the matter has no meridional motions, the only motion 
of the matter is a circular one that is represented by 
the angular velocity, the angular velocity is constant as 
seen by a distant observer at rest and the matter can be 
described as a perfect fluid. To study the rotating stars 
the following metric is used 

ds 2 = -e^ +p) dt 2 + e 2a {dr 2 + r 2 d9 2 ) 

_l_ e (7-p) r 2 s j n 2 Q(d(j) — udt) 2 (32) 

where the gravitational potentials 7 , p , a and u are func¬ 
tions of polar coordinates r and 9 only. The Einstein’s 
field equations for the three potentials 7 , p and a can be 
solved using the Green’s-function technique [l3j and the 
fourth potential ui can be determined from other poten¬ 
tials. All the physical quantities may then be determined 
from these potentials. Rotational frequency of stars is 
limited by Kepler’s frequency which is the mass shed¬ 
ding limit. For very compact stars such as neutron stars 
the Kepler’s frequency is very high and can go up to mil¬ 
lisecond order [lj, [l5| whereas white dwarfs being about 
thousand times bigger in size and much less dense, Ke¬ 
pler’s frequency is very small and one may safely use 
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the zero frequency limit fl6| to the Einstein’s field equa¬ 
tions. Obviously, at the zero frequency limit correspond¬ 
ing to the static solutions of the Einstein’s field equations 
for spheres of fluid, the present formalism yields the re¬ 
sults for the solution of the Tolman-Oppenheimer-Volkoff 
(TOV) equation jl7, |lg| given by 

dP{r ) G [e(?') + P(r)][m{r)c 2 -f 47rr 3 P(r)] 
dr ~ C 4 r 2[l 2 Gm(r) ] 

L 7*C 2 J 

where e(r) = p(r)c 2 and m(r)c 2 = / £(r')d 3 r' 

Jo 

which can be easily solved numerically using Runge- 
Kutta method for masses and radii. The quantities e(r) 
and P(r) are the energy density and pressure at a radial 
distance r from the centre, and are given by the equa¬ 
tion of state. The mass of the star contained within a 
distance r is given by m(r). The size of the star is de¬ 
termined by the boundary condition P{r) = 0 and the 
total mass M of the star integrated up to the surface R is 
given by M = m(R ) [lj|. The single integration constant 
needed to solve the TOV equation is P c , the pressure at 
the center of the star calculated at a given central density 
Pc- 


V. RESULTS AND DISCUSSION 

Recently, there are some important calculations for 
masses and radii of magnetized white dwarfs using non- 
relativistic Lane-Emden equation assuming a constant 
magnetic field throughout which provided masses up to 
2 .3-2.6 Mq [20 |, a mass significantly greater than the 
Chandrasekhar limit. However, because of the struc¬ 
ture of the Lane-Emden equation, pressure arising due 
to constant magnetic field do not contribute while for 
the general relativistic TOV equation case is not the 
same. Moreover, the EoS needed to be fitted to a poly¬ 
tropic form. In order to derive a mass limit for mag¬ 
netized white dwarfs (similar to the mass limit of ~1.4 
Mq obtained by Chandrasekhar (2lJ for non-magnetic 
white dwarfs), the same authors, under certain approx¬ 
imations, have been able to reduce the EoS to a poly¬ 
tropic form with index 1 + 1/n = 2 for which analytic 
solution of Lane-Emden equation exists (0(£) = sin£/£ 
where p = p c 6 n with p and p c being density and cen¬ 
tral density, respectively) and avoiding the energy den- 
sity £b = and pressure Pb = ^£b arising due to 
magnetic field by assuming it to be constant throughout, 
they were able to set a mass limit of 2.58 Mq [22], f23| . 
In the present work, we have calculated masses and radii 
of white dwarfs by solving the general relativistic TOV 
equation both for non-magnetic and magnetized white 
dwarfs using the exact EoS without resorting to fit it to 
a polytropic form. 


A. Chandrasekhar limit for White Dwarfs 

We verify Chandrasekhar limit (2lJ for masses of 
white dwarfs by actually solving TOV equation for non¬ 
magnetic white dwarfs. The masses and radii of such 
white dwarfs are listed in Table-I. It is interesting to note 
that considering a very high central density of 3.343x 10 10 
gms/cc for / = 1 white dwarfs, one can asymptotically 
reach the Chandrasekhar mass limit. It is important to 
mention that beyond this density at ~ 4.3xlO n gms/cc, 
the neutron drip point [24], the nuclei become so neutron 
rich that with increasing density the continuum neutron 
states begin to be filled, and the lattice of neutron-rich 
nuclei becomes permeated by a sea of neutrons. In Table- 
1, masses and radii of non-magnetic white dwarfs as a 
function of central density are provided. In Fig.-l plot 
for masses of non-magnetic white dwarfs is shown as a 
function of central density whereas in Fig.-2 mass-radius 
relationship of non-magnetic white dwarfs is provided. 
These results for non-magnetic white dwarfs do conform 
to the traditional Chandrasekhar mass-limit. 

TABLE I: Variations of masses and radii of non-magnetic 
white dwarfs with central number density of electrons which 
can be expressed in units of 2xl0 9 gms/cc for mass density 
by multiplying with 1.6717305X 10 6 . 


n e (r=0) 

Radius 

Mass 

fnT 3 

Kms 

Mq 

l.OxlCT 5 

917.87 

1.3904 

5.0X1CT 6 

1126.83 

1.3905 

4.0X1CT 6 

1202.53 

1.3896 

3.8X1CT 6 

1220.55 

1.3893 

3.6X1CT 6 

1239.80 

1.3890 

3.4X1CT 6 

1260.43 

1.3887 

3.2X1CT 6 

1282.65 

1.3883 

3.0X1CT 6 

1306.67 

1.3878 

2.8X1CT 6 

1332.78 

1.3873 

2.6X1CT 6 

1361.33 

1.3866 

2.4X1CT 6 

1392.75 

1.3859 

2.2X1CT 6 

1427.62 

1.3850 

2.0X1CT 6 

1466.67 

1.3839 

1.8X1CT 6 

1510.90 

1.3825 

1.6X1CT 6 

1561.69 

1.3809 

1.4X1CT 6 

1621.01 

1.3788 

1.2X1CT 6 

1691.86 

1.3761 

l.OxlCT 6 

1779.00 

1.3724 

8.0xlCT 7 

1890.72 

1.3673 

6.0X1CT 7 

2043.29 

1.3594 

4.0xlCT 7 

2275.36 

1.3457 

2.0xlCT 7 

2721.16 

1.3138 

l.OxlCT 7 

3233.63 

1.2692 

l.OxlO -8 

5482.58 

1.0051 

l.OxlO -9 

8721.75 

0.5949 
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FIG. 1: Plot for masses of non-magnetic white dwarfs as FIG. 2: Plot for mass-radius relationship of non-magnetic 

function of central density. white dwarfs. 


B. Super-Chandrasekhar White Dwarfs 

As mentioned in the beginning of this section that u 
like non-relativistic Lane-Emden equation, pressure aris¬ 
ing due to constant magnetic field does contribute to 
the general relativistic TOV equation. Presence of high 
constant magnetic field do not provide valid solutions 
to the TOV equations. Hence, we present stable so¬ 
lutions of magnetostatic equilibrium models for super- 
Chandrasekhar white dwarfs with varying magnetic field 
profiles which is maximum at the centre and goes to 10 9 
gauss at the surface of the star. This has been obtained 
by self-consistently including the effects of the magnetic 
pressure gradient and total magnetic density in a general 
relativistic framework. Nevertheless, we have also per¬ 
formed calculations corresponding to very high (single 
Landau level) and high (multi Landau levels) magnetic 
field which is constant throughout the star in order to 
compare with the results from solutions of Lane-Emden 
equation described above, but for these cases we have to 
ignore the explicit contributions from energy density £b 
and pressure Pb arising due to magnetic field. Results 
of such calculations are provided in Tables-II & III for 
magnetized white dwarfs with single and multiple Lan¬ 
dau levels, respectively. 

Now we perform the actual calculations with varying 
magnetic field including the effects of energy density and 
pressure arising due to magnetic field in a general rela¬ 
tivistic framework. The variation of magnetic field [25j 
inside white dwarf is taken to be of the form 

B d = B s + B 0 [l - exp{-^(n e /n 0 ) 7 }] (34) 

where Bd (in units of B c ) is the magnetic field at elec- 


TABLE II: Variations of masses and radii of uniformly mag¬ 
netized white dwarfs with central number density of electrons 
which can be expressed in units of 2xl0 9 gms/cc for mass 
density by multiplying with 1.6717305X 10 6 . The minimum 
magnetic field Bdmin corresponding to the central density re¬ 
quired to make single Landau level throughout is listed in 
units of B c . 


n e (r=0) 

Radius 

Mass 

Bdmin 

fm -3 

Kms 

M 0 

in units of B c 

5.0xl(T 6 

592.28 

2.4521 

253 

4.0xl0 -6 

636.54 

2.4508 

218 

3.0xl0 -6 

698.11 

2.4461 

180 

2.0xl0~ 6 

792.71 

2.4204 

138 

l.OxlO -6 

989.84 

2.4149 

86.5 


tronic number density n e , B s (in units of B c ) is the sur¬ 
face magnetic field and no is taken as n e (r= 0)/10 and 
/3, 7 are constants. Once central magnetic field is fixed, 
Bo can be determined from above equation. We choose 
constants /3 = 0.8 and 7 = 0.9, rather arbitrarily by us¬ 
ing unequal non-unity values, which provides stable so¬ 
lutions of magnetostatic equilibrium models for super- 
Chandrasekhar white dwarfs. Nevertheless, the mag¬ 
netic field is not taken completely in ad hoc manner, 
because central and surface magnetic fields once fixed 
the variations of its profile do not alter the gross results. 
Moreover, we have kept maximum central magnetic field 
strength at 10 B c which is 4.414 x 10 14 gauss, near to the 
lower of the maximum limit suggested by N. Charnel et 
al. [26} and surface magnetic field ~ 10 9 gauss estimated 
by observations. In Table-IV the results of these realis¬ 
tic calculations are listed. In Figs.-3,4 plots for masses 
and radii of magnetized white dwarfs are shown as func- 

















2 4 6 8 10 

B dc (in units of B 0 ) 

FIG. 3: Plot for masses of magnetized white dwarfs as a func¬ 
tion of central magnetic field. 

TABLE III: Variations of masses and radii of uniformly mag¬ 
netized white dwarfs with central number density of electrons 
which can be expressed in units of 2xl0 9 gms/cc for mass 
density by multiplying with 1.6717305X 10 6 . The magnetic 
field Bd (< Bdmin for the central density) is also listed which 
generates multiple Landau levels. 


n e (r=0) 

Radius 

Mass 

B d 

fm“ 3 

Kms 

M© 

in units of B c 

4.6736x10“® 

1149.77 

1.3925 

1.5 

3.5147x10“® 

663.58 

2.4491 

200 


tions of central magnetic field. Present calculations esti¬ 
mate that the maximum stable mass of magnetized white 
dwarfs could be ~3 M & . These results are quite useful in 
explaining the peculiar, overluminous type la supernovae 
that do not conform to the traditional Chandrasekhar 
mass-limit. 


VI. SUMMARY AND CONCLUSION 

In summary, we have considered a relativistic, degener¬ 
ate electron gas at zero temperature under the influence 
of a density dependent magnetic field. Since the electrons 
are Landau quantized, the density of states gets modified 
due to the presence of the magnetic held. This, in turn, 
modifies the EoS of the white dwarf matter. The pres¬ 
ence of magnetic held also gives rise to magnetic energy 
density and pressure which is added to those due to de¬ 
generate matter. We hnd that the masses of such white 
dwarfs increase with the magnitude of the central mag¬ 
netic held. Hence we obtain a conclusive result that it 
is possible to have electron-degenerate magnetized white 



2 4 6 8 10 

B dc (in units of B c ) 

FIG. 4: Plot for radii of magnetized white dwarfs as a function 
of central magnetic field. 


dwarfs, with masses signihcantly greater than the Chan¬ 
drasekhar limit in the range of M 0 , provided it has 
an appropriate magnetic held prohle with high magni¬ 
tude at the centre as well as high central density. 

TABLE IV: Variations of masses and radii of magnetized 
white dwarfs with central number density of electrons which 
can be expressed in units of 2xl0 9 gms/cc for mass density 
by multiplying with 1.6717305X 10 6 . The maximum magnetic 
field Bdc at the centre is listed in units of B c whereas the sur¬ 
face magnetic field B s is taken to be 10 9 gauss. 


n e (r=0) 

Radius 

Mass 

Bdc 

fm“ 3 

Kms 

Mq 

in units of B c 

4.674017x10“® 

1285.91 

1.4146 

1.5 

4.673846x10“® 

1344.46 

1.4236 

1.75 

4.674209x10“® 

1349.45 

1.4339 

2.0 

4.675374x10“® 

1388.04 

1.4906 

3.0 

4.672188x10“® 

1438.94 

1.5731 

4.0 

4.670830x10“® 

1503.64 

1.6863 

5.0 

4.678118x10“® 

1581.27 

1.8353 

6.0 

4.677677x10“® 

1663.86 

2.0217 

7.0 

4.665741x10“® 

1758.40 

2.2601 

8.0 

4.661657x10“® 

1954.44 

2.8997 

10. 


To date there are about ~250 magnetized white dwarfs 
with well determined helds jj] and over ~600 if ob¬ 
jects with no or uncertain held determination [ 23 , j28| 
are also included. Surveys such as the SDSS, HQS and 
the Cape Survey have discovered these magnetized white 
dwarfs. The complete samples show that the held dis¬ 
tribution of magnetized white dwarfs is in the range 
10 3 -10 9 gauss which basically provides the surface mag- 





























netic fields. However, the central magnetic field strength, 
which is presumably unobserved by the above observa¬ 
tions, could be several orders of magnitude higher than 
the surface field. In fact, it is the central magnetic 
field which is crucial for super-Chandrasekhar magne¬ 
tized white dwarfs. However, the softening of the EoS 


accompanying the onset of electron captures and pyc- 
nonuclear reactions in the core of these stars can lead to 
local instabilities which set an upper limit to the mag¬ 
netic field strength at the center of the star, ranging from 
10 14 -10 16 gauss depending on the core [26j composition. 
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